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INTRODUCTION 

This progress report presents results of research carried 
out. under Grant- NGL-34r-002-047 during the period February 1 to 
May 31, 1970. 

As a part of the efforts to use standardized symbols, and 
notations in research under the grant, a list of notations related 
to wire grids is presented. 

The report contains next a paper which appeared- in the 
Transactions of the IEEE on Microwave Theory and Techniques in 
the January issue of 1970, and which deals with the H-guide with 
laminated dielectric. Lamination is one of the approaches to 
reduce the dielectric losses in the H— guide- which contribute a 
major part of the attenuation of- this guide. The study shows 
that a reduction of the ^ loss tangent by a factor of about 10 can 
be achieved by this approach. The improvement can be obtained 
by using thin strips of low-loss material with high permittivity 
separated by layers of air or by layers, of a, foamed low-loss 
dielectric. 

In an- effort to reduce, the wall' losses in the waveguides 
under investigation, such as H— guides, fence -guides, and, re- 
flector guides- by proper treatment of the • conductor surfaces, 
the effects of the surface roughness and surface structure on 
the- attenuation are being investigated. Optical methods are used 
to measure appropriate parameters for the surface roughness. 
Various degrees. of surface roughness are obtained by polishing 
the surfaces on papers of various, grain sizes. Once calibrated, 
the optical setup permits simple but reliable measurement of the 
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surface roughness. The waveguides under investigation, ' in turn-, 
permit accurate determination . of the effects of surface roughness 
on the Q-value and on, the surface resistance of the conducting 
surfaces. The measurements are' expected to permit determination 
of optimum characteristics- obtained by proper surface treatment. 

A detailed study of the error sources of Q-value measure- 
ments- was- carried out which indicated the possibility of obtain- 
ing- improvements by refining the measurement procedures and by 
eliminating frequency jitter resulting from vibrations entering 
the microwave oscillator used for the frequency measurement. 

A modulation- oscillator for the X-band, source was then added 
to the measurement setup which permitted simultaneous determina- 
tion -of the half-power frequencies. By the, use of water cooling 
and. vibration- free support of the X-band circuitry, frequency 
jitter of the X-band source was reduced considerably. As a con- 
sequence of these .measures , the .-measurement error was reduced 
to about +0.2%. The changed- procedures- and the improved setup 
are described- in part 4 of the report. 

Discrepancies were observed at the comparison of measured 
and computed, values of the attenuation of the fence guide. This 
made a review of the relationships used for the computation - of 
the- characteristics of wire grids- necessary . It was found that 
the relationships derived and described in. the literature are no 
longer valid for dense wire grids as. used in the fence guide. New 
relationships were., derived for- the surface impedance which adequately 
describe the properties of wire grids with dimensions useful for the 
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fence guide; Equations for lossless and, lossy wires were found. 
Derivations and the relationships are described- in part 5 of this 
report. 

Part 6 presents the x~esults of measurements of the characteris- 
tics of the fence guide carried out on a, prototype section of • guide 
with a length of- 175 mm. Two such sections with different spacings 
of the wires were fabricated and one of them investigated during 
the report period. Prior to the measurements, the measurement 
setup and the- field probes, were tested and improved to give in- 
creased accuracy and improved reproducibility-. Longitudinal 
and transverse field distributions were then measured and the 
results evaluated. The -results are ■ presented in part 6 of the re- 
port. Reduction. of the diameter of the compensated probe gave 
considerably improved results- in comparison with those obtained 
in previous measurements. The- measurements continue on the 
second section of fence guide with increased wire spacing. 

An experimental study was carried out. on a shorted section 
of- open reflector guide previously investigated theoretically, 
(Progress .Reports of June 15 and October 15, 1969). The- shorted 
section represents a resonator whose Q-value- and- internal- field 
distribution were measured. It was found that improved resonance 
conditions and improved Q-values (as high as 30,000) can- be ob- 
tained by reducing the spacing of the reflector surfaces. To 
study the effect of a variation of the. distance on the Q-value 
and the field distribution, one of the - reflector surfaces was 
made, movable. This permits determination of the characteristics 
of the resonator under nonfocal conditions,. Preliminary results 



of the . measurements are presented- in -part 1 7 of- the .report . 

As a part of the study, of- dielectric loading of H-guide and 
♦ 

fence guides, the electromagnetic field-matching method has been 
studied. It was- applied for finding the characteristics of a • 
parallel-wall, waveguide with a centrally located metallic- bar. 

In comparison with other approaches, the method gave satisfac- 
tory results. Computer programs. were developed for the use of 
the method and. for the comparisons. The- study is described in 
part 8 of the report. 
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STANDARDIZED NOTATIONS RELATED 
TO WIRE GRIDS 

. Impedance of wire grid itself without losses 
(reactive part) . 

. Contribution to wire-grid impedance by con- 
duction losses of wires.' 

. Total impedance of wire grid. 

. Surface impedance taking into account radiation 
through the grid. 

. Spacing between . axes • of wires-. 

. Diameter of wires. 

Geometry parameter. 

Normalized reactance. 

Normalized resistance. • 


[Z w /2 o. = 3P + q] 
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H Guide with. Laminated Dielectric Slab 


FREDERICK J. TISCHER, fellow, ieee 


Abstract — A consideration of Uic field distribution and the charac- 
teristics of an H guide with a laminated dielectric liar -composed of an 
arbitrary number of dielectric strips separated by air layers is presented. 
Low-loss wave modes with the electric field intensity predominantly paral- 
lel to the sidewalls arc analyzed. Approximate thin-layer equations arc 
derived which indicate methods for improvements of the attenuation of 
the H guide. 

Introduction 

T HE H guide is a hybrid waveguide in which the fields 
are concentrated in one transverse direction by sur- 
face-wave guidance and are confined in tile other direc- 
tion by parallel conducting planes. The cross section has the 
fofin of an H with the conducting side walls forming its 

Manuscript received January 10, 1969; revised July 7, 1969. This 
work was supported in part by NASA Grant NGR 34-002-47. 

The author is with the Department of Electrical Engineering, North 
Carolina State University, Raleigh, N. C. 27607. 


vertical legs and the dielectric slab the horizontal' bar. The 
guide has low conduction losses in the walls for wave modes 
with the electric field strength predominantly parallel to the 
sidewalls, similar to those of the waveguide with circular 
cross section for TEio waves. Since the major contributions 
to the attenuation of the waveguide result from the dielectric 
losses, lamination of the dielectric bar was proposed [l] to 
reduce these dielectric losses. Computation and measurement 
of the characteristics of a double-strip //guide confirmed the 
feasibility of improving the loss characteristics by this 
approach [2], [3], [6]. 

This paper presents the derivation of the field distribution 
and of characteristic data for an H guide with the dielectric 
bar subdivided into an arbitrary number of dielectric strips 
separated by air layers. The resulting relationships can be 
drastically simplified in the case of thin strips. Under these 
conditions, the slab structure may be considered as -an arti- 
ficial dielectric with an equivalent permittivity e Tcn and loss 
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Fig. I. H guide with laminated dielectric, (a) Cross section of the guide, 
(b) Laminated-slab structure. 


langcnt tan 5 oq . Diagrams are presented which show the 
attainable improvements in the loss characteristics of the 
juidc resulting from the lamination. Low-loss wave modes 
ivith the E field in the center of the guide parallel to the side 
.vails are considered. 


strength, which is parallel to the boundary surfaces (TM), 
and the otiicr field components have in customary notation 
the following form for waves propagating in the positive 
2 direction: 


II v w = { yl,. exp [-y Xn (x — x A -i)] 

+ B„ exp [+Yxn(a; — t„_i)] j c~ lk ' z , 
, x j dII v M 

Ex M = J! — = (k/we)///*), 

oit dz 


(la) 

(lb) 



we 


dll./ ni 

dx 


= 0y*»/“0 { An exp [— Yx„(.c — .r„_i)] 
— B n exp [+ 7 ,„(® - x»_i) ] } 


(1c) 


The propagation constants are related in air by 


Yxn cc x ] ct x ~ — &o 2 j Ic x > h oj ko — w\/ eo^o, (2a) 
and in dielectric by 


Method of Derivation 

The method of deriving relationships for the laminated 
H guide is based on the concept that the fields in a parallel- 
.vall waveguide can be considered as those of superimposed 
jlanc waves bounced back and forth between the conducting 
, vails. Applying this approach to the present case, the analy- 
sis can be carried out in two steps. First, surface-wave propa- 
gation along an infinite slab of a laminated dielectric is con- 
sidered, and then two plane surface waves are superimposed 
o form the field within the //guide. 

The structure under investigation is illustrated in Fig. 1(a) 
or the example of a slab composed of six dielectric strips. 
Since the distributions of the field' components are symmetri- 
cal and antisymmctrical with regard to the horizontal center 
>Iancs, it is sufficient to consider the fields in the upper half 
jf the guide, as if a perfectly conducting electric or magnetic 
vail separated the upper from the lower half. 

The outlined method of derivation hence suggests as a 
irst step the derivation of relationships for the surface-wave 
iropagation along the laminated infinite dielectric slab, 
ts structure, representing only the upper half, is indicated 
n Fig. 1(b). It consists of N layers consisting alternately of 
lielectric and air. In reality, the slab may be composed of 
ligh-permittivity dielectric strips separated by layers of 
bam material with low permittivity and low loss. The dielec- 
ric strips have equal thicknesses d„ characteristics e r , and 
eparation's d a . 

tVave Propagation Along a Laminated Dielectric Slab 

Wave propagation along a laminated dielectric slab was 
:onsidcred by Mathis [4] who used a transmission-line ap- 
>roach for the determination of the propagation charactcris- 
ics. This approach will be used in a modified form for 
leriving the field distribution in the stab structure of Fig. 
(b). The wave fields arc described in rectangular coordinates 
is solutions of the wave equation. The magnetic field 


T-Tfi jk XJ h x ~ — 1c r A, 'C A" q \/ (2b) 

Summation of the equations for a x 2 and k x 2 , (2a) and (2b), 
yields 

+ K 2 = A’o 2 (e r - 1), (3) 

which is one of the basic equations describing the propaga- 
tion characteristics of the surface waves. 

The field distribution in region 1 («=1) for magnetically 
symmetric wave modes is found from (1) by introducing the 
boundary condition E C =Q for x=0. In customary notation, 

/fj, fl) = 2A\ cosh a z x, 

E x m = (AyW 0 )2.'li cosh a*x, (4) 

EA l) = — j(a x /ue 0 )2Ai sinh a x x. 

The corresponding equations for electrically symmetric (odd) 
modes (H u —0 for x— 0) are found to have cosh a z x and sinh 
a x x of (4) interchanged. The fields in dielectric region 2 are • 
given by 

H v m — exp [— jk x (x — Si)] 

+ B 2 exp [+jfc,(a: - *0]}, 

E x ™ = fa/ue^R,™, • ’ . (5) 

E t m = — (fcx/co6 0 «r) { A 2 exp [— jk z (x — Xi)] 

- B -2 exp [+jk x (x - xj)]}. 

The constants A 2 and B« can be found in terms of/L for 
even and odd modes by equating the tangential field com- 
ponents //„ (1> , E z w and H„ {2) , ES 2) at the boundary at x—xi. 
Equating the field components at the next higher order 
boundary and subsequent repetition yields the constants A n 
and B n and after insertion into (1) the field components in 
terms of the amplitudes of the field components at the 
lower boundary in the center of the guide at.*=»0. In the 
uppermost air region which extends to -v— «, the field com- 
ponents are, since. j5, v = 0, given by 
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U v m - Ax exp [-a* (a; - z.v-i)], 

E * ( ' v > = (fc./coe* (6) 
— j(.a x /wt 0 )Atf exp [~a z {x — - x,v-i)J. 

The derived equations show that the field distributions are 
sinusoidal in the .v direction in the dielectric strips and are 
described by hyperbolic functions in the air layers. The hy- 
perbolic functions indicate the surface-wave character in the 
air layers. In the air region above the laminated structure, 
the field decreases exponentially toward infinity, assuming 
infinitely wide strips. 

The multiple-layer structure according to Fig. 1(b) and 
tiie form of (1) suggests derivation of recursion equations 
which express the field amplitudes at each upper boundary 
of a layer in terms of those at the next lower boundary. Con- 
sidering only components in the directions transverse to x, 
we find from (1) 


Use of the elements of the matrix T facilitates expressing the 
transverse wave impedance at the upper boundary of the 
structure in terms of the fields and impedances at x=0: 

AV' V) 2WV»(0) + T^EA'KO) 

Zj? (12) 

H.w TuH» w (Q) + TuEA»(0) 

Substitution of the impedance value Z iV from (6) and of the 
impedance values at x-0 from (4) yields the characteristic 
equations for the laminated structure: 

cx T' 

(E, = 0 at x = 0) : Z* < — j = — — , (13a) 

weo Tn' 

for even modes, and • 

a* 2% 

(H„ - 0 at x = 0) : Z N ~ j — = - — -, (13b) 

C06o F 12 


• = //„ tn ~ 1) (.T n _ 1 ) cos 

1 

+ 7 — E/ ,, - ll (x n - 1 ) sin k z d n , 

E z M (x„) — ZpJI^^fxn-^j sin M, 

+ E t (n ~ l) (x n - 1 ) cos k z d n . 


(7a) 

(7b) 


The equations are similar to those for current and voltage 
at the input and output of sections of transmission lines. 
Hence, the field distribution in the x direction can be con- 
sidered by transmission-line methods as indicated by Mathis. 
According to the transmission-line concept, a transverse 
impedance can be defined as Z (n) =E r M /HF n) which can be 
expressed by familiar relationships as a function of Z (n-1) 
with Z p „ k z , and d„ as parameters. The latter describe the 
characteristics of the nth layer, where 


7 _ 1 7 . 7 _ P<>. 

' J V « — - "Of ^0 — ‘ * 

€ r Kq 60 


Similar equations involving hyperbolic functions govern the 
relationship between the various quantities in the air layers. 

It is convenient to introduce matrix representation for 
the consideration of the multiple-layer structure. Writing 
(7) in matrix form, a transmission matrix [r] relates the 
field quantities. For the nth dielectric layer, it is given by 


for odd modes. Inspection of (9), (11), and (13) shows that 
the characteristic equations describe relationships between 
k x and a z for even and odd modes. In combination with (3), 
they facilitate the determination of k z and a z . The transcen- 
dental equations (13) can then be written in a general form 
as 

kx fc,o(.&X) €r, dtj da) ( 1 4) 

where d, and d a are the thicknesses of the dielectric strips 
and the air layers, respectively. The exact evaluation of the 
characteristic equations can be carried out for large numbers 
of layers by computer. An impedance-matching process by 
the use of the Smith chart outlined by Mathis [4] can be 
used if the number of layers is small. Approximate relations 
derived in the next section are satisfactory for many practical 
cases. 

Relationships for Thin Dielectric Slabs 

The possibility of simplification of the derived equations 
for thin layers is explored next since it can be expected that 
the slabs of the laminated H guide satisfy this condition. 
Assumption of thin layers, when a z d a , allows the 

following approximations: 

cos k z d t , cosh a z d a 1 ; sin k z d t ~ k x d t ; 


v 

cos k z d n j — sin 


[r<»>] =1 

L jZ p , sin k-xdn cos k z d n 


Transmission in all the layers below can be taken into con- 
sideration by writing 


sinh aj a c* a z d a . 

The approximations are introduced into (7) and (9) and- 
into two representative transmission matrices T M and T {n ~ l) 
for a pair of layers, one containing air and the other con- 
taining a dielectric medium. By matrix multiplication, a new 
matrix can now be found that describes the transmission 


r Hu<n) ] = rv'tmirrtn— ol . rvcol [^V^ 0 )"] through a uniform single layer equivalent to the combined 

J ‘ ^ J ’ ’ ' L-* L^* a) ( 0 )J" transmission through the two layers. The new transforma- • 

tion matrix is hence given by 

Taking into account the total number of layers between- ATiv_i 


and x=0, a cumulative transmission coefficient [T] can be rytm] _ • 
defined for the Inmlmued structure: 1 J L 


'T tv^-i 

J 


\i U I 1 "] “ • • • [Tm]. (11) 

_ [• 

TiS"- 

-y u (n-i> 


Li 21 1 22- 

L Tu M 
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Substitution of the approximate values of the matrix ele- 
ments into (9) and (15) and carrying out the matrix multipli- 
cation yields 


T i,«» as 1 ^ 1 - — , (16a) 

€r 

Tn (D) ~ SreqC d, + d a ) ~j — (ed, + d a ), (16b) 
Zo Zo 

€roq 


. Zo/K* \ 

- j " ( df ocx'd a ) j 
ko\ (r / 


Tti (D) ^ 1 ^ 1 -f a T -d t d c 


(16c) 

(16d) 


where the permittivity, propagation constant, and thickness 
of the equivalent single layer are denoted by c rc<l , k Icn , and 
d t +d a , respectively. Evaluating (16b) yields an approximate 
expression for the equivalent permittivity, 


1 + p/e, 

^ e, — — » 

1 + V 


(17) 


y 

with p=d a /d, denoting the air-to-dielectric thickness ratio. 

Since the characteristics of the equivalent layers are equal, 
the total laminated dielectric structure can be replaced by a 
uniform dielectric slab with equal total thickness and having 
a relative permittivity « rcq . For thin dielectric layers, the 
laminated dielectric H guide can thus be treated as a guide 
with a single slab with an equivalent permittivity given 
by (17). 


Thin-Strip Loss Tangent 

Caution is required at the determination of the loss tan- 
gent of the equivalent uniform dielectric representing the 
thin-strip structure. We observe that e r „i of (17) actually is 
the real part of the complex equivalent permittivity. It was 
obtained by considering wave propagation and the field dis- 
tribution in the x direction and involved E. and H v [see (7)]. 
The imaginary part €r«,"(e rc q =€,<*,'— j^') which describes 
tire dielectric losses must be found by a more detailed con- 
sideration of the losses and attenuation of the equivalent 
layer, since the major contribution to the dielectric losses 
results from Ex which did not appear in the equations from 
which e rC q' was derived. This component (£ x ) has its maxi- 
mum value in the center of the slab at x=0 for low-loss even 
modes. At this location the component E- is zero. Due to the 
sin 4 distribution of E z> its contribution to the losses can be 
neglected in the central region where the dielectric slab is 
located, arid only E x taken into consideration. 

Determination of e rc<1 " will be based on the concept that 
the attenuation of the equivalent single layer is equal to 
.that caused by dielectric losses in the combination of the 
original dielectric and air layers. The compared values of 
attenuation arc which represent the ratios of dis- 

sipated and transmitted powers per unit width of the layers. 
The transmitted power is derived from the z component of 
the Poynting vector, 


(18a) 


(18b) 


Substituting H v for E z from (5) and o-,=weoe/', we find, for 
the two layers with the dielectric giving a contribution only, 

^-d t , ' (19a) 

to«o (e,)" 

and for the equivalent single layer, 

Pdh ; r- j \’h v I* — (c u + do).' (i9b)' 

«60 («roq E 

Equating the unprimed and primed ratios of pa^/pu^m yields 

d, W' 

d, -f- e/da 

and subsequently for the equivalent loss tangent 

tan S 

tan - e K ” / e r «/ ^ • (20) 

1 + pe/ ■ 

It is interesting that the validity of (20) can be checked by 
considering a capacitor with a laminated dielectric with the 
boundary planes parallel to the capacitor plates. The equa- 
tion for the equivalent loss tangent is found to be identical 
to that of (20) for low-loss materials. 

Fields in the Laminated H Guide 

The fields in the H guide are found by superposition of 
two plane surface waves propagating along the infinite lami- 
nated slab discussed in the preceding sections. One of these 
waves travels toward the sidewall incident at a finite angle 
from the normal ; the other, the reflected wave, travels toward 
the center of the guide. The geometries of the directions of 
waves are shown in Fig. 2. After a change of notation, the 
field components previously denoted by H v and E x are now- 
called the transverse components H t and E, of the plane 
surface waves. The component E, becomes the longitudinal 
component 2s { . Hence, 

II v — > II p, E x —* Ei ) E - — > Ei] fcj — > ki. (21) 

By use of the geometric relationships indicated in Fig. 2, we 
express the field components in the coordinate system of 
the H guide [see Fig. 1(a)] as a superposition of incident 
(e+^y) and reflected waves (e~^»»); hence, omitting the 
common term er> k ‘ c , the field components are 


S. = i- Re [EJI ?* ], 
which gives for the double layer 

Ptran, = * | | 2 ~ ^0 + A 

*o W / 

and for the equivalent single layer 

k Z 

Ptran/ = £ | II „ | 2 -J (d, + d a ). 

&0 €roq 

The dissipated power is derived from 

Pdiss “ 2" { 
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Fig. 2. Geometric relations of the field components of incident and 
reflected waves at the wall of the guide. 


H z — H 4- tr iL * v ) = 2 HtQcJki) cos k v y, 

II v = II i{k.-/k,)(eJ^v — = 2II t (k t /ki)] j sin kjj, 

E- = Ei{c ll * u — e~ }k » v ) — 2 Eij sin k v y, 

E v = - Bi(k u /ki)(e *#> + e~^«) ‘ ' (22) 

= — 2Ei(k u /ki) cos k u y, 

E t = Ei(k,/ki) (c> l » u — cr’ 1 *”) — 2Ei(k./kj)j sin k„y. 

Substitution of the field components for the infinite lami- 
nated slab described by (4) through (7) into (22), after taking 
into account the changes of notations expressed by (21), 
. gives the field distributions in the corresponding layer within 
the H guide. For the nth layer which contains dielectric 
' material, the field=components are 

II, M {x) = y \U n c(x) + jVnS(x)] COS kylJCT > k ‘ z , 
k[ 

■ kr 

Hy M (x) = j — [U n c{x) + jF n s(.r)] sin kyijer**, 
lit 

[Er M (x) - ■ — — [U n c(x) +jF„s(a;)] sin k v ye . 

OHatr \ .. ( 23 ) 

E u M (x) •= — — y [F„c(a:) — jU n s(x)] cos k u yc~ j7!,z , ■ - 
W6o hi 

Er M {x) = ~j~l~ y [F n c(a;) ~ jU n s(x)\ 
cvenfr ki 

■ sin k,jye~ jkt ‘. 

Tile functions c(.v) and s(.v) stand for cos kj(x—x„-i), find 
. sin kr(x— .v„_i), respectively. The constants U„—(A n +B,i)/ 2 
and V „ = (/!„ — B„)/2 arc proportional to the field components 
Er and E, at the lower boundary of the region Tire cor- 
responding equations for an air layer are found by substitut- 
ing hyperbolic functions for the trigonometric functions c(x) 
. and j(.v), < r = 1, and —ja x for k x . 


Mode Analysis 

Confinement of the waveguide fields in they direction- by 
the conducting walls is associated with an infinite spectrum 
of higher order wave modes. The boundary conditions, E t — 0 
at y=0, a, require that k vm —mir/a. Since the longitudinal 
propagation constant ki is related to its components by 
ki i =k u i +k z t (Fig. 2), the propagation constants in the vat- 


ious directions are interrelated by 


kifer — k x ~ + k um ' -f- k zm - 

(24a) 

in dielectric, and 


ko~ = - a * 2 + kyj + hJ 

(24b) 


in air. In a waveguide of proper width, only the fundamental 
mode is propagating and the higher order modes are evanes- 
cent with the fields decreasing exponentially in the z direction 
with negligible phase shift. 

The constants a s and k x which describe the field distribu- 
tion and the surface-wave character in the X direction can 
be computed from (3) and (14): 

A* 2 = V(<, -1) - ax 2 , ' (25a) 

kr* = fcAax, €r, d a , d t ), ■ (25b) 

which are the characteristic equations for the infinite lami- 
nated slab. The relationships are not affected by-the confine- 
ment by the two parallel walls and are valid for the H guide 
also. Excitation of the even and odd modes with different 
distributions in the x direction depends on the operational 
frequency, the permittivity of the strips, and the geometry 
of the laminated structure. The fundamental magnetically 
odd wave mode (even with regard to Ex) can always be 
excited since its existence is frequency-independent. The 
higher order modes with multiple periodicities -in the ,v 
direction in the slab and all magnetically even modes have 
cutoff frequencies as lower limits for their existence. The 
propagation characteristics are basically similar to those of 
surface waves along an infinite solid dielectric slab [5]'. The 
similarity becomes, particularly pronounced in the case of 
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thin layers of the dielectric and air for which the approxi- 
mate relationships derived previously can be used. Denoting 
the single or multiple- solutions- of the equation pair (25) by 
A' r „ which have to be found for the specific operational 
conditions, the propagation constant in the z direction be- 
comes 







( 26 ) 


Equation (25b) has to be solved for thick layers by computer ; 
for a small number of strips, it can be solved by the imped- 
ance-matching process outlined by Mathis. In the case of 
thin layers the substitution by a uniform dielectric slab with 
an equivalent permittivity can be applied. Basically the 
excitation of waveguide modes in the laminated waveguide is 
similar to that in the H guide with a single solid slab de- 
scribed previously. 


Improvement of Loss Characteristics by Lamination 



Fig. 3. Equivalent permittivity of the laminated dielectric slab. 


The study of the two-strip H guide has shown that the 
attenuation has a general trend to decrease witii decreasing 
thickness of the dielectric slabs [2], [3 j. The same trend can 
be observed in the case of the laminated guide with thin 
dielectric strips. It can hence be expected that a considerable 
reduction of the attenuation can be achieved by lamination 
of the dielectric slab and by reducing the thickness of the 
dielectric strips. Due to these trends the thin-strip approxi- 
mation should give satisfactory results in a consideration of 
the attenuation of the laminated H guide. 

The major contribution to the attenuation of an H guide 
at millimeter waves of nonexcessive size results at present 
from dielectric losses. Information with regard to the reduc- 
tion of tile attenuation can be gained by evaluating and com- 
paring the dielectric losses of various structures disregarding 
the effects of differences of the geometries on the wall losses. 
If the equivalent permittivities of the compared slabs are 
kept constant, the disregard of geometric changes is well 
justified. The relationships derived for thin layers will be 
used for such an evaluation next. 

The dielectric losses in the guide tire directly proportional 
to the loss tangent (tan 5) of the medium of the strips. A reduc- 
tion of tan 5 is hence the primary source of a reduction of 
the attenuation. The equation which shows the feasibility 
for such a reduction is (20), 

tan 5 n , = tan 5/(1 -f- e r p) ; p = d a /d t . 

It should be noted that this equation involves primarily the 
electric field component E x which contributes to the trans- 
mission of power along the guide. It has a maximum ampli- • 
tude in the center of the guide where the dielectric is located. 
The other components E u and E. are zero at x=0' [see (4)] 
and can thus be disregarded- in the approximation. The 
second important equation involves the equivalent permit- 
tivity . •- ; 

— e r (l + p/«r)/( 1 ’ + P).' ‘ 

Graphical representations of the two relationships are shown 



4. -Reduction factor of loss tangent of the laminated dielectric 

slab ( 7 = tan Soo/tanS.u.p). 

in Figs. 3 and 4. The preceding considerations indicate that ' 
a considerable reduction' of the attenuation of the H guide ; 
can be expected at millimeter waves if tan 5 cq can be reduced ' 
by choice of proper materials and geometry p of the lami- 
nated structure while keeping constant. The equations- 
indicate that dielectrics with high permittivity at low dielec-r 
trie losses andTamination with a high value of p are best- 
suited for, a reduction of attenuation. The comparison of loss 
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TABLE I 



Single-Slab 

Laminated-Slab 

Dielectric 

Parameters 


Parameters 



tan S 

P 


tan 

Styrofoam 

1.03 

1 .5X10-* 




Rexohtc. 

2.53 

7X10-<* 

7.5 

1.16 

3.5X10- 6 

Beryllium Oxide 

6.6 

4.4X10-* 

3 

2.4- 

2.6X10- 5 

(99.5%) 

Aluminum Oxide 

9.9 

2.5X10" 5 

5 

2.5 

5X10-’ 

(99.9%) 

Eccoceram 

90 

1.0X10“’ 

10 

i 

1 9 

1.1X10“* 


* At 24 GHz, others at 10 GHz. 


tangents and permittivities in Table I shows that equivalent 
loss tangents in the guide obtained by lamination are con- 
siderably superior to those of dielectrics in the single-slab 
guide. It seems that a reduction of the loss tangent by one 
order of magnitude is possible. 

Conclusions 

Wave propagation in an infinite laminated slab has been 
•considered as a first step of an analysis of the laminated 
H guide. Equations for the field distributions in the various 
regions of a multiple-layer structure have been presented and 
characteristic equations for the total structure derived. 
Matrix representation was found most useful for describing 
the characteristics of the multiple structure for even and odd 
wave modes. The infinite-slab consideration represents the 
first part. In the second part, the H guide fields are derived 
by superposition of two infinite-slab surface waves incident 


and reflected from one of the two parallel conducting side- 
walls. This approach simplifies considerably the derivation 
of the field equations. The subsequent analysis of the wave 
modes shows that the wave modes have basically the same 
characteristics as those of the single-slab guide. 

Thin-layer approximations allow considerable simplifica- 
tions of the derived equations. They facilitate simple con- ■ 
sideration of the laminated structure by an equivalent single 
layer of a medium with an equivalent permittivity and loss 
tangent. Emphasis is then placed on an evaluation of the 
equivalent loss tangent in a thin-layer structure since it indi- " 
cates clearly the methods by which the attenuation of the 
guide can be reduced particularly in the millimeter-wave 
region. The equations, the graphs, and a comparison of the ' 
equivalent loss tangents for various materials show that use - 
of low-loss high-pcrmittivity materials and combination 
with thin-strip lamination give considerable reductions of the 
dielectric losses which in turn lead to considerable improve- 
ments of the attenuation of the H guide. 
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RESULTS OP LIGHT REFLECTANCE MEASUREMENTS 
FROM ROUGH SURFACES 

Introduction 

In a study of the effects of surface roughness on the loss 
characteristics of waveguides and cavities at millimeter waves, 

optical methods are used for the determination of the surface 
roughness. 

It is well known that light is reflected from metallic 
surfaces but the reflection phenomena are still being investi- 
gated. If the metallic surface is polished to a mirror-like 
finish, light is reflected at an angle to the normal equal to 
the incident angle. However, for a surface containing irregular 
ities, the reflected light contains both specular and diffused 
parts. The specular light is reflected as previously stated 
but the diffused light is scattered- in all directions. As -the 
roughness of the surface increases, then for the same incident 
light the intensity of specular reflection decreases while the 
diffused increases. For near normal incident light it has been 
observed that the diffused light intensity with angular dis- 
placement about the specular angle appears to have a Gaussian 
distribution. From experimental data the spread or variance 
of this distribution is observed to be dependent on the "rough- 
ness" of the surface. This paper is a description of the pre- 
paration of the surfaces, the equipment used to measure the 
reflectance, and preliminary data reduction. A microscopic 
analysis of the surface is in progress and will be presented 
later. 



Sample Preparation 


14 


Since the final goal - is to determine the contribution- of 
surface roughness. to attenuation- in waveguides, -the -metal chosen 
was copper. A sheet. of 1/8" copper was cut into lj by 2 inch 
sections. This size was selected for ease of handling as- well 
as constraints imposed by the reflectance- measuring equipment. 

Grinders manufactured by Buehler Ltd. - were available and 
could be used to prepare. the samples. There. were- two types 
of grinders; one consisted of a set of abrasive papers mounted 
on hard, flat surfaces and' the other having the same type of 
mounting but a continuous flow of water immediately washed 
away the loose abrasive material and fragments of the metal. 
These metallurgical grinders, both the dry and wet, were used 
for preliminary preparation in polishing of- metal. surfaces. 

The copper surfaces were prepared in the following manner. 
First, the sample was moved forth and back on the coarsest grade 
of paper until grooves- were .present in only one. direction. It 
was important that the grooves appear. in only one direction 
since the finer paper can only remove the grooves created by 
the next coarser paper. The force on the samples was held to 
a minimum to prevent crushing the abrasive particles which 
could cause large variation in- the grooves and cause smoother 
surfaces than intended. The- force was only that to keep the 
samples in a position to create' parallel grooves. ‘ The sample 
was then -rotated ninety degrees and placed on the next. finer 
paper. The process was continued until the desired degree of 
polishing was obtained. In our investigation, six grades of 
abrasive paper were used. The grain size of each abrasive paper 
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is given in table- 1. 

After each grinding the samples were rinsed under running 
water to remove loose material from the surface. A stream of 
air with sufficient force was blown across the surface to 
rapidly remove residual water and thus minimize water spots. 

The surfaces were carefully handled to prevent finger prints or 
other stains and were stored in a cabinet- until reflectance 
studies could be made. 


Reflectance. Setup 

The reflectance measurement setup is illustrated in Fig. la,b. 

The light, source is .a Spectra-Physdcs Model 124 Helium-Neon Laser 

o 

capable of supplying fifteen milliwatts of power at 6328 A. An 
adjustable attenuator with degree indicators, along the circum- 
ference is used to regulate the light intensity that strikes 
the test sample. The relative attenuation can .be determined 
by the angular displacement from a reference when applied to 
the sin square law. A continuous monitor of the source inten- 
sity is accomplished- by a beam splitter which deflects part of 
the light, into a Spectra-Physics . 401 C- power meter. If varia- 
tions. should occur during an experiment the attenuation can be 
adjusted to correct the- change. A flat, circular disc with a 
mounting stand located in. the center is used to secure the test 
sample. The- five inch radius disc has degrees marked off around 
the perimeter and- with aid of markers, the angles of incidence 
and reflectance can- be measured. An extension arfoi which pivots-, 
about the- center of the disc holds a light detector which con- 
sists of an RCA 7102 multiplier phototube and' associated circuitry. 



Power 

meter 



Fig. la Reflectance measurement setup 
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Fig. lb Photograph of reflectance measurement setup 
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Table 1 

Abrasive Paper Particle Size 


Grade 


Grain size in microns 


Emery 2 62 
Emery 1 38 
Grit 240 48 
Grit 320 30 
Grit 400 24 
Grit 600 15 
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The photo tube is enclosed in a light tight cover with a one- 
eighth by one inch slit. The extension arm maintains the slit 
at a constant distant from the center of the disc and at a 
constant height above the plane of the disc. The power for the 
multiplier is supplied by a Hewlett-Packard Harrison 6516A DC 
power supply and the output of the detector is displayed on a 
Doric Integrating Microvoltmeter. The assembly, with the ex- 
ception of the power supplies and the voltmeter, is mounted on 
an optical rail which. permits- accurate alignment. of the various, 
components. The entire setup is located in a .dark room whose 
walls have been sprayed with a non-glossy black paint. 

Three sets of experimental data were taken from the setup. 
For the first- experiment, the angle- of incidence and incident 
light were held constant, and in a plane perpendicular to sur- 
face, the reflected light was measured at angles from the sur- 
face jiormal between 0 and 90 degrees. During the second experi- 
ment, the incident angle was held constant and the light re- 
flectance was measured at the specular angle for various inten- 
sities of incident light. For. the final experiment the incident 
light was- held constant and, the incident angle was varied between 
10 and 90 degrees while the reflectance was measured at the 
specular angle. These. measurements were. taken for each of the 
test samples and is given in- the next section i 

Results 

The data were plotted as three separate graphs and are shown 
iri Figs. 2, 3, and 4. Figure 2 is - a plot of the normalized re- 
flectance versus angular displacement, from the surface normal for 



Normalized reflectance 


1.0 . 



10 20 30 40- 50 60 ’ 70 


Angle of .reflectance in degrees 

Fig. 2. Reflectance of rough surfaces averaged and normalized 
with regard to maximum versus angle of reflectance 
for constant angle- of incidence 



Light detector output in volts 


12.0' I 


600 



Power input in milliwatts 

Fig. 3 Light intensity (relative) reflected from rough surfaces 
versus incident light power. Angles of incidence and 
reflectance were equal and kept. constant. 


f— 1 


Reflectance in, volts 





Normalized reflectance 



10 - 


Angle of. reflectance in degrees 

Fig. 5 Individual measurement results which form the 
basis ^ of the averaged reflectance curves of 
Fig 
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the case where the incident light and angle are -held constant. 
The curves for each test sample are the results - -of averaging of 
three locations on the sample. The* locations are - the center 
and 1/4 inch on each side of the center. Figure 5 shows the 
variation and. how the different locations, relate to- the average. 
The curves appear symetrical about the specular angle and re- 
semble Gaussian curves.. The coarser the abrasive paper the 
broader is the spread or variance- of the curve. 

Figure 3 shows a- plot of specular reflectance versus in- 
cidence light for a constant angle, of incidence. .Each set of 
data. gives a linear curve. For these cases- the coarser the 
abrasive paper the smaller the amplitude of the reflected light 
and the greater the slope. 

Figure 4 illustrates the dependence of roughness on the 
angle of incidence. The graphs are plots, of the specular re- 
flection versus the incident angle for the case of constant inci- 
dent, light. The expected result would be that at some angle the 
surface would reflect .more, light than hearer the .normal and would 
increase its reflectance until the grazing angle was obtained. 
This was verified because- the finer the abrasive paper used the 
sooner the increase began and the greater the slope of the curve. 



Discussion: . 


The model, used in this section is shown. in- Fig. 6. 



Fig. 6‘ 

The surface lies in the xy plane and the incident light is in 
the yz plane and makes an angle ij; with surface normal. Note 
that the reflected light does not reflect entirely in the yz 
plane. From this a relationship for the reflected light is 
formulated: 

R == I M(a o ,e/U) G(a>a,X ,if>, 6, 4>) 

where R is the intensity of the reflected light- (reflectance) 
I is the intensity of incident light 
M- is a variable- dependent on- the material 
properties (cr Q ,e,y) 

G is variable- dependent on the geometrical, properties 
■T}ve,4> are angles as given, in Fig. - 6. 
a, a are constants for the description of the 
statistical properties of the surface 
X is the wavelength of the incident light 
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Since the same material, copper, is used in all cases, the 
quantity M is a constant factor. For the' first -set-of data, Fig. 
2, the- variables are the angle- e for each sample. The height 
variation is associated with various grades of -abrasive paper. 

It is noted that for each test sample- the resulting graph re- 
sembles a Gaussian curve but the spread is different and depends 
on the grain size of the abrasive paper. It is reasonable to 
assume that the larger grain size would produce greater surface 
roughness. The data confirm this assumption. 

The next set of data shown in Fig. 3; has the one. obvious 
conclusion that the polishing with the finer grain size gives 
the greatest reflectance. The slopes of the curves confirm 
this. Using the proposed general relation for reflectance, the 
slope is given as 


- M(cr,-e,u )G(a,a,X ,\p, 0 , <j>) 

The results indicate the slope- is smaller for the rougher surface. 
From this and, the previous experiment the variation of the • 
function is . assumed to have a Gaussian form with the spread in- 
creasing- with -increasing surface roughness. 

H-. Davies considered light, reflection by a rough surface 
with a Gaussian density function for surface .variation and cal- 


culated the reflectance for two cases; a slightly rough and a 
very rough surface.. The two results are: 


Slightly rough: 


R - I K. 


a 2 a 2 


(cose + cosijj) 4 exp{ - 


2 2 
7T a 


I(sin0cos<j> - sinij)) 


+ sin 2 0sin 2 <J>] } 
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Very rough: 

R = I K v ~ ( cos 0 + cos 4>) 2 exp{~ sL (sin9cos'(i) - sin^) 2 *sing 6sin?<j> j 

° 2 2 a 2 ( cos0 + cost ] >) z 

where K g and K^. are constants and- a' and - a 'are -constants de- 
scribing the statistical properties of the ■ surfaces . The 
equations indicate the general trend of the experimental data. 

The various parameters are needed to relate the equations to 
the experimental results. 

The final graph illustrates the dependence of the surface 
roughness on the angle- of incidence. The Rayleigh criterion 
for surface is given, by assuming light strikes -a maximum and 
minimum height of the surface and is reflected in the same 
direction (see Fig. 7). The phase difference between the 
two paths is. 

A a . 47Th , 

A<p = — cos ip 

where h is the total height, difference and ip is the angle of 
incidence. If A<p — tt> then cancelation occurs and no energy 
flows, in this, 'direction. If Acp — 0, it reflects specularly 
and is considered smooth. A value between the two extremes 
is more often used as- indication of- roughness, then for a 
smooth surface 

h < r . 

8 cos 4* 

From the simple formulation • it is obvious that as the angle 4> 
approaches ninety degrees moderately rough surface appears 
smooth. This effect is easily observed in Fig. 4. The 
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rougher surfaces require large incident - angles before they 
begin to show characteristics of a smooth surface. 



Appendix 

Rayleigh Criterion t 
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2 



consider two rays incident on a surface with height variation 

of h and at an incident angle -The path difference between 

the two rays is 

£r = 2h cos ip 

«no therefore the phase difference is 

^ 6r * cos i^i 

For^ small phase' difference- the rays -are in, phase -and. reinforce 
hut xf the- phase difference is tv they cancel. - Since this means 
no energy flow in this direction, it - is- concluded that it has 
been redistributed in other directions or the -surface is con- 
side red rough. Then- between the two extremes .is a transition 
region; id « */2 has- been chosen as veil as other -values as an 
indicator of roughness. Using this value then for a smooth 
surface 

h < .. . . 1 

8 cos, p ' 

From this a surface tends. to be smooth for either hA + D or 
tfi V' 10°, 
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AN IMPROVED Q-VALUE MEASUREMENT TECHNIQUE 
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Introduction 

In- the June 15, 1969 progress report a description of a 
frequency measurement setup was given (1). The method was used 
in Q-value measurements to measure half-power frequencies and 
resonant frequency by- observing beat frequencies with a single- 
frequency oscillator. The system proved to give Q-values of a 
resonator to a typical accuracy of ±150 parts out of 8000. A 
more, stringent requirement was created by- a surface-resistance 
study where the order of the maximum Q-value- change was ±50. A 
method using a modulated oscillator for the frequency indication 
which gives this accuracy is described. The improvement. gives 
Q-values of ±50 parts out of 8000. An up-to-date description 
of the setup is also. given. 

Circuitry of the General Test Bench- 

The bench setup as shown. in Fig.-l has a 35 GHz A-band 
Klystron as the signal, source. Power is supplied to the A-band 
signal source by a Narda 62A1 power supply. The Klystron signal 
is fed through an- isolator, an adjustable attenuator and a lOdb 
directional coupler. From the main guide of the coupler,- the 
signal travels through an- isolator to a precision attenuator and 
then into the resonator. under investigation. The secondary guide 
of the lOdb directional- coupler is - fed into a 3db directional 
coupler which serves as a power divider. The-, signal of the 
secondary arm of the 3db coupler is fed through an adjustable 
attenuator and a cavity-type frequency meter to .a- tunable detector 
mount. The main arm of the 3db cdupler is connected to a waveguide 




Waveguide attenuator 

mixer 


Fig. 1 Frequency measurement setup 


CO 
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switch. The signal from, one of the switch positions is fed 
through an adjustable attenuator into a waveguide mixer. The 
other switch position may be used for other experiments (presently 
not in use) . 

Frequency Measurement 

Since no frequency measurement equipment is presently 
av silable for frequencies up to 100 GHz, a special setup was 
developed for these measurements . 

For swept frequency operation the repeller plate of the 
A— band klystron is modulated with . a 60 Hz sawtooth waveform. 

This makes the klystron signal vary in frequency proportionally 
to the amplitude of the sawtooth, voltage. The sawtooth voltage 
is also used as the input to the horizontal deflection plate of 
the HP 1200B oscilloscope. This allows the horizontal deflection 
of the oscilloscope to be as an approximation proportional to 
frequency. When the klystron signal is displayed on the oscillo- 
scope (as a function of the repeller voltage) , each position of 
the trace represents a frequency.. In order to know the frequency 
which corresponds to each horizontal position of the deflected 
beam the heterodyne principle is being used. 

A c.w. signal of known frequency (within the range of our 
display, which we can find by use of a cavity- type frequency 
meter) is mixed with the A-band signal. The A-band signal is 
changing linearly in frequency if we keep the amplitude of the 
modulating sawtooth relatively small. The- known c.w. signal .is 
the result of frequency multiplication within the mixer and 
originates in an x-band source that can be tuned to have its 
harmonic frequency within the range of that of the swept A-band 
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signal. The x band signal is fed into the frequency multiplier 
mixer. The x-band frequency is continuously monitored on a high 
accuracy frequency counter. 

The harmonic of the x-band signal forms a beat- frequency 
signal with the A-band signal. If the interference signal is 
superimposed on the mode pattern. observed on the screen of the 
oscilloscope, a zero beat pattern can be seen. This pattern 
actually is - a sine wave with a frequence equal to the difference 
between the instantaneous A-band frequency- which is swept linearly 
by the sawtooth and the harmonic of the x-band signal frequency. 

A zero beat occurs at a specific position of the beam when the 
known harmonic of the x-band signal is equal to the unknown 
A-band signal.. By this, the horizontal-, deflection of the oscillo- 
scope {sweep proportional to the frequency- of the A-band signal) 
can be calibrated in terms of the frequency with counter accuracy. 

Circuitry for , Frequency Measurements 
A Model 54 OA function generator by Exact Electronics is used 
for generation of a, sawtooth, wave . with a frequency of approximated 
60 Hz. This signal is- applied to the horizontal input of the HP 
1200B oscilloscope, and -also to. the external modulation jack of 
the 62A1 klystron power supply for the A-band klystron. This 
makes the frequency of the signal of the A-band varying in the 
form of a sawtooth linearly with time. To keep. the operation 
linear the sawtooth has to be. relatively small. By applying this 
same waveform to the horizontal deflection plates of the oscillo- 
scope, the display is then a function of frequency. On channel B 
of the oscilloscope the mode pattern of the A-band klystron can 
be monitored. With a cavity-type. frequency meter the approximate. 
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center frequency of the klystron mode, can be- found. This approxi- 
mate frequency- divided by the number of a convenient harmonic 
yields the frequency of the x-band. source. - ' 

The c.w. -x-band. signal- is generated by a Varian-X-13- x-band 
klystron- and is fed into a Demornay Bonardi BBGD-35G crystal 
multiplier-mixer. The frequency of the x-band source is monitored 
by a HP 52 45A. frequency counter. The other input signal to the 
crystal- mixer is the sawtooth -modulated A-band- signal. Variable 
shorts on the crystal mixer permit adjustments for optimum fre- 
quency multiplication and. mixing. The crystal current, is ' not to 
exceed 50 jia.' 

The output signal of the crystal mixer is amplified in an- 
AC amplifier (lOHz-to 1.5 MHz) and fed into an RC-fi-lter. The 
output. is capacitively coupled to the line- connecting the crystal 
detector for monitoring the klystron mode to channel B of the 
oscilloscope. On the screen of the oscilloscope, the mode 
pattern is then displayed, with ■ the interference signal superim- 
posed ©nr it. The control knob for the reflector voltage on the ' 
X-13 signal- source .power supply (HP 716B Klystron Power Supply) 
can be used to bring the two frequencies closer together. ' The 
zero beat between. the two spikes of the interference pattern 
corresponds to the condition that both frequencies - are equal. 

The x-band source is constantly being monitored by the accurate 
HP 5245A frequency counter. Multiplication of- the monitored fre- 
quency- by the number of the harmonic gives the unknown frequency 
at the horizontal position of the -beam indicated by the zero beat. 

The measurement of a, segment of the scope trace' corresponds 
to a difference in frequency. If the x-band. source reflector 
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voltage is modulated (HP Signal Generator 606A)*. The> indica- 
tion on the scope has the form as shown in. Fig. 2b . Because 
we are using the fourth harmonic, one fourth of the frequency 
of the modulating signal is the difference frequency on. the scope 
of the first sideband from the carrier. ' The beat-frequency 
spectrum displayed on the- oscilloscope is used to measure the 
half-power frequencies of the resonance curve. 

Q-Value Measurements 

The resonator to be measured is inserted into the experi- 
mental branch of the general setup. The signal from the A-band 
klystron is coupled into the resonator and the output resonance 
curve is displayed on channel B of the Hp 1200 B oscilloscope. 

On channel- A the mode, pattern of the klystron is displayed. By 
adjusting the positioning of the traces. and the attenuation of 
the signals, one signal is placed on the top of the other so 
that the peak of resonant curve touches the klystron mode curve 
(see. fig. 2a ). -The resonant curve is. increased by removing 
3db of attenuation making the intersection of the two curves 
be at the half-power points of the resonant- curve. Than taking 
the- beat-frequency spectrum- of the- indicator when the reflector 
voltage. of the x-band klystron is modulated and aligning the- 
first sidebands with- the intersections, the half -power- fre- 
quency difference is found.. We now have a display that appears, 
like. that in Fig. 2b. Since. the carrier- frequency as displayed 

*The x-band. signal- consists of the carrier and side-frequencies 
AF apart, where AF.is the modulation frequency-. 



'ig. 2 Dual beam displays on oscilloscope for a Q-value measurement, 

(a) Initial setting. 

(b) Resonance input increased by 3 dB. 

A. Cavity resonance curve. 

B. Klystron— mode . curve . 

B . Klystron-mode curve "with superimposed interference pattern. 

C. Zero level. 
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Fig. 3 Monitor console for frequency measurement setup 
a X-band power supply 
b Dual trace oscilloscope 
c Function generator (sawtooth wave) 
d Frequency counter 
e Microwave circuitry 
f Modulation signal source 
g X-band Klystron 
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on the oscilloscope is. that of the x-band source,: the correct 
frequency (because we use the fourth harmonic) is • four times 
that of the x-band signal. The first sideband frequency- must 
be doubled to get. the total frequency difference (note that 
the first sideband frequency display- on the scope. is one-fourth 
the frequency of the modulating signal) . With points F 1 and F 2 
corresponding to the half -power frequencies, the loaded Q is. 

V'W- 


References 

1. R. A. Kraft and. F. J. Tischer: . Frequency Measurement 
Setup for Millimeter Waves, . Progress . Report- - Grant 
NGR 34-002-047/52, . pp 28-38 (June 15, 1969). 



CHARACTERISTICS OF DENSE WIRE GRIDS 


Abstract 

As a supplement to preceding considerations , this report 
deals' with- the characteristics of wire grids- with narrowly spaced 
wire's-. The presently used equations are reviewed from the viewr 
■point of their application- at decreasing spacing between the 
wires. Modified equations are developed for lossless - and lossy 
grids for the -region of spacing between zero - and values where 
customary equations' can be applied. 

Introduction 
- L 2 

In .preceding '-reports ' the characteristics- of wire grids 
were, considered to find relationships to be. used for the 'study 
,and' the design of the - fence guide; At .the derivations, equations 
were used from- a. few major articles and, the results com- 

pared with those. obtained- by experiments. Satisfactory agreement 
was obtained. At the continued study of the fence guide, less 
spacing between. the wires of the grids -was applied recently. This 
had the consequence that the known equations are no longer appli- 
cable and that a modification- of the previously used relationships 
became necessary to- allow computations in ' .the . region of, decreasing 
spacing. Derivation of the modified equations and the resulting 
relationships for lossless and lossy wire grids are' presented in 
this report. 

Review of. Relationships ' 

Two basic relationships derived by various authors are known 
and. can be used for the computation of. the surface impedance of 
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wire grids. One of these was derived directly for grids with 
round wires and has the following form 

- 1 Z — - 5 - . (!) 

The other equation represents a modified 'relationship derived - 
originally for grids- composed of strips. At the modification^ 
the surface of the circular conductors was assumed to equal 
that of the strips. The- resulting equation is 

= J ^ csc IP * (2) 

We observe that substituting the argument l/(ird/s) for the 
cosecant, which is valid for large spacings,’ .makes the two 
equations become identical. Since in the- present case. the. 
spacing and the wire diameter usually are much smaller than the 
wavelength, a correction factor, originally introduced by Mac- 

3 

Farlane , is disregarded. 

Considering the application of these equations to the fence 
guide, it was found that the second equation gives a larger range 
since it can be- used for smaller spacing than- Eq. (1) . The lat- 
ter [Eq. (1) ] gives a negative impedance if the term s/ird becomes 

smaller than unity.- It was also found that Eq. (2) gives better 

2 

agreement . with the experimental results . 

The equations can be rewritten in- somewhat different form 
by combining the geometry parameters to form separate- terms We 
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write hence 


Z P /2 0 

z p /z 0 



(3a) 

( 3b ) 


It becomes- advantageous to substitute- for the geometry term s/d 
a- new variable t = s/d - 1. This variable actually indicates 
the width of the gap between the .wires normalized with regard 
to the wire diameter d. The new variable becomes- zero when 
s/d approaches one and the wires touch each other. 

It is interesting to inspect plots, of the second factors 

of Eqs.. (3) which are shown, in Fig. .1- as functions' of t. Curve. 

A represents- the second, factor of-Eq. (3a) and B the corfespond- 

ing term, of Eq. (3b).- The graph also contains the values of 

this factor obtained by the evaluation of, measurements of the 

surface. impedance- in the X-Band regions described in previous 

2 

reports. They are indicated by. D and E. We observe that the 
two curves A. and B- approach the abscissa at finite- values of 
t. This- indicates that the corresponding relationships are in- 
correct, for small spacing and' not- valid since the impedance can- 
not approach the value _ zero . while a- gap exists -between the 
wires.’ A modification- of the relationships - is hence necessary 
to satisfy the condition that the impedance approaches gradually 
and. monotonously zerd. for t. = 0. 

As the next. step, an empirical method is chosen to find- an. 
approximate relationship to represent the- second factor in the 
equation for the wire grid impedance in terms - of t for small 



Relative values of geometry, factor 



L 2 3 4 5 


Geometry parameter, . t = ■ s/d- - 1 
Fig. 1- Geometry factors for dense wire grids. 

A. (s/d) An (s/ird) . 

B. { s/d) An esc ( ird/ s ) . 

C. . {s/d - 1) 2 /4; t < 4. 



44 


values of this quantity. The' conditions for the choice of a 
useful relationship are’ the following:' The relationship should 
be simple, it should give zero for p = 0 and should gradually 
become Eq. (2 ) for increasing values of- p. A relationship 
based on t was chosen to accomplish this. The^ relationship 
is shown, graphically as. curve C in- Fig. 1 with a = 2 as- a 
satisfactory, choice-. It is given • mathematically by 



L | Q 


We observe that the quadratic relationship valid for s/d < 5 
(t < 4) assumes the value of. the second factor in Eq. (3b) for 
s/d =• 5 or t = 4; The measured value indicated by E is now- 
very close to the theoretical ‘ value- obtained by the new relation- 
ship. Since- the value- Of the -second factor of Eq. (3b) is 2.67 

- t 

for t = 4 (s/d = 5).,. the- new approximate equation for the im- 
pedance for small • spacing becomes 

Z p /Z 0 = j (f ) O.I67t 2 (£) 0Ab7 (f- if- (4) 

Effect of- Conduction- Losses, 

It can be assumed. that the customary relationships for 
taking into consideration the losses- on the- wires are.no longer 
valid at small spacing. These relationships take into account 
the losses by assuming that the current along each wire is uni- 
formly distributed around -the wire and’ by. adding to. the surface 
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impedance of the. grid a term 5 

where R g is the surface resistance of the grid material. 

There are two -aspects which have, to be., taken- into considera- 
tion for small spacing. First, the distribution in- circumferen- 
tial direction of the- current along, the wires. is -no longer uni- 
form and is such that- the major part of the current flows on the 
input half of the circumference. This- is- the' case' since, due 
to the. considerable reflections -,by the dense grid, the fields 
and, the surface current are low on- the - output- s^ide,. Hence, only 
one-half of the circumference ird [in the denominator of E<j. 
is practically available to carry the- grid- current. Furrth'ernu^re , 
the current will not be uniformly distributed on. the input half 
of the circumference of the grid wires. It will be a- maximum- 
in the center where the magnetic field' lines are dense and low 
on the sides of the cylindrical wires where the surfaces of 
adjacent wires face each other. . The local ■ contributions to the 
losses, by the parts of the surfaces on' the sides will be very 
small ( since these contributions are proportional- to the square 
of the relative local surface currents and local ' magnetic fields. 
A satisfactory approximation is offered by assuming a sinusoidal 
distribution of the current around the input half of the wire 
circumference. This allows, simple determination- of the maximum 

current density in the center. The- integration of the. current 

t cur rent: 

over the input half of the circumference gives the to tali' per wire 
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Fig. 2 Circumferential- distribution of 
longitudinal current on wires. 
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resulting from the incident waves. We can write 

- T — 

'"'max ~ ^ 


0 *) 


since the average . value of a, sine wave is 2 /tt.. Based on this 
current. distribution, we can .subsequently determing an equivalent 
surface- resistance R' 1 according to 

rd 


D oy Rfi Jl '4 = f ^ x) M X j 


0 

where. the distribution of J 2 (x) is given -by the dotted line in 
Fig. (2) . Evaluation of the integral yields 


n 2 r, 7 ^ iL - ~] 9 ‘ 1 IT d 2 j[f! 

D av % ^ Z max 2 ^ Z = ( Z ' 4 


It follows 


T 


R ' - 
s 8 


R s = 1.23 R 4 


a) 


Taking into- consideration- that only one-half of the circumference 
of the wire is available to carry this current, we find for the 
contribution to the impedance of the wire grid 


2 


w<r 


TTd/Z 


0+j) 


2.46 




This indicates- that for a, closed wire grid the effect of the 
losses increases by about a factor 2.46 in comparison with wide 
grids . 

The variation of the- constant, 2 „ 46 . with- increasing spacing 
can.be taken into account by substituting for it a term which 
linearly decreases. with increasing t. . The assumption- is made 
thatj for a spacing given by t = 4^ the constant has decreased to 
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a value 1- giving the- equation- the customary form expressed by 
Eq. (4&). The relationship which satisfies these conditions is 


given by 

- 0,365 1 ~ 


It follows that the approximate equation for the- contribution of 
the conduction losses to the surface impedance has- the form 


R. 


« ( 2Ab - t) 


Of) 


w<r ' T d 

The equation is- valid- in the. region '0 < p < 4 (1 < s/d < 5). 

The- total approximate .surface impedance of the wire grid 

in this- region is hence-. 


°- i67 2 o 4- (i- | ) 2 +c 2 - £,6 -°' 36Fi )||fi + -j)' / 

ax act 


w-J -oj -\ - 'ttc V”'*'- (Q) 

The derived relationship will be used. until mo re ^equations will 


become, available. It can be expected, however, that mo’re 
accurate equations will not drastically change the relationships. 
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FENCE GUIDE MEASUREMENTS 


Abstract 

Measurements of the field distribution in 175 mm long 
fence- guide sections- are discussed and an evaluation of -the 


parameters is made. 
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In the previous progress -report, (NGL- -.34-002-047 of February 
15, 1970) electromagnetic fields measurements _ were., carried out 
at 35 GHz on a, fence- guide model 47.5 mm long. The measurement, 
setup was- also included. At present measurements are carried 
out on a longer, section of fence guide. 

Two 1-75- mm long fence- guide- sections were fabricated using 
copper wires- of diameter d w = 0.862 mm, spacings- s w = 1.17 mm 
and 1.Q16 mm on a dielectric sheet witha thickness.2d = -0.80-mm 
and e r = 2.53 (Rexolite) . The prototype sections are shown in 
Fig. 1. The measurements were made at 35.0 GHz and an improved 
capacitive monopole with compensation was- used. as ■ a field, probe. 

A thin copper sheet was- used as a, short- in some cases of the 
measurements and the output response' was . recorded with the help 
of an X-Y recorder. 

The VSWR pattern with load end shorted (Fig'. 2) was measured 
over- the, full, 175 mm- length of the. fence ' -guide. The-. minima 
over the entire length were -precisely located and a mean separar 
tion of 5. 442. mm found. "The guide wavelength for- the . fence guide 
at. 35 GHz computed- from- the average distance between the minima 
is - 10.884 mm. 

g 

The field decay was- measured in the x-direction- away- from 
the dielectric, with the fence guide short, circuited, at the 
locations, of first, . second' and the. fourth maximum and. up to- the 
height of the fence wires. The -exponential • decay of the electric 
field was- plotted . on a linearized- db scale..- (Fig. 3). If averaging 
technique, is used, the lines are straight and, fairly parallel to 
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Fig. 1 Measurement setup for the study of an open 
cylindrical-reflector resonator. 






Relative amplitude in.dB 



Fig. 3 Exponential decrease' of field strength 
inside the fence- guide 
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each. other. The slope of the lines gives the exponential decay 
factor a varying somewhere in the range of 1.67-dB/mm to 1.80 
dB/mm. The- corresponding decay factor as measured on the Standard 
H-guide has an. average value. of 1,53-dB/mm. 

Field measurements were made to compare' the-. level of the 
field. .strength, inside {Fig. 4) at the center and outside (about 
3 mm away from the fence) with a short circuited load. The field 
strength, is indicative of the leakage -of field through the- wires 
of. the fence. The- maximum field measured outside was at least 
36 dB below- the maximum- field inside. 

The variation of the x-component of the electric field in 
the transverse direction- was- measured at the location of the 
first three maxima from the short circuit for various- heights 
above, the dielectric sheet (Fig. 5) . The variation has a sine 
form as- was expected. It showed some irregularity in- the form 
for the third half-period. 

Figure. 6 shows' the field distribution measured in transverse 
direction at the . level- of the upper edge of the. fence wires. The 
parameter is the distance-, above the . upper edge. 

In an effort, to obtain, some. idea, of the attenuation- in the- 
longitudinal direction, VSWR measurements- were' carried out using 
a tunable -matched load.- The pattern, obtained- had a reasonable- 
trend for the envelope and, the attenuation was- preliminarily 
computed at a standing-wave- ratio 1.05 for the load'. The value 
so obtained, was .066 dB/mm. The measurements- are being continued. 
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STUDY, 0F OPEN RESONATORS - WITH CYLINDRICAL REFLECTORS 

Abstract 

An experimental investigation of a non-confocal open 
resonator, with, cylindrical reflectors is reported. The 
resonator was operated at 36,8GHz. Measurements were made of 
the Q-value- and of the. field distributions ' as . a function of re- 
flector spacing. Mode numbers .were . obtained from these plots 
and an attempt was- made to find relationships • between the mode 
numbers m,n-,p, the spacing * b 1 between the reflectors, and the 
Q- value. 


Introduction 
1 2 

In- several- papers ' related to the study of- open resona- 
tors for millimeter waves, confocal systems, (distance of separa- 
tion between the , reflectors equal to the radius - of curvature of 
the reflectors) which are considered optimum from the point of 
view of diffraction losses were investigated, . The- possibilities 
of instabilities- of this- system were reported, 3/4 

Investigations made in this laboratory on two confocal 
resonators, with- cylindrical reflectors, manufactured to normal- 
machining tolerances- gave unsatisfactory operation and relatively 
low Q-values,. However, regions of considerably higher ’Q* of the 
order of 3.0,000 were observed with- the- same- reflectors for re- 
duced distance between the reflectors (r/2 < b < r) Based on 
these results it. can be concluded that deviations, from confocal 
geometry may be preferable. The tolerances - at the fabrication 
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of the two reflectors and also the uniformity of- spacing between 
them may not have to be as rigid, as required for- the confocal. 
case. A detailed study of the properties of . non-confocal systems 
is in progress . 


Experimental 'Study 

The- main • purpose of the. present investigation is to study 
the . properties of the open resonators with- cylindrical reflecters 
over a wide range of spacing between. them and to find- out- whether 
there exists any relationship between the Q-value, the geometry 
parameters, and the mode numbers. Field distributions are 
measured inside- the cavity at resonance, the mode. numbers deter- 
mined, and the. Q-values found. The setup used for the investiga- 
tions, is .shown in Fig, 1. The geometry of the cavity is indicated 
in • Fig. 2 . . 


Results 

The distance of separation- was- varied within the range 

0.5 < r/b <1. A manual,. system was used to move one of the re- 
flectors with respect to the other in- such a way that the aper- 
tures of, the reflectors remained parallel and- their symmetry planes 
coincide; In two ranges the observed Q-values were of order of 
30,000. The- distance between. the reflectors in these. two ranges 
were 

1 . b ~ r/2 + 0 » 5 ■ cm 

2 . - b =* r - 0". 5 cm 

A considerable mode interaction and. irregular resonance patterns 
were observed for other spacings. The resonance pattern of some. 
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Fig. 1 Prototype sections (L = 175 mm) of fence guide 




X,Y, z.: Rectangular Coordinates Axes 

m: Mode number along ' x'axis 

n: Mo de> number along- axis 

p; Mode number along ' z - ' axis 
Q: Quality factor 

bs Spacing between reflectors 
r: Radius- of curvature of the reflectors 

Fig. 2. Geometry of resonator 
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of the high- ' Q' - modes wer.e not symmetric, and no correlation 
between the Q-values and spacing was- found. 

Measurement of Field Distributions. 

The aperture distribution as well as the transverse dis- 
tribution were obtained by moving a compensated probe • inside the 
cavity and. using the output of the probe- to drive a plotter. 

The plots, are shown ; .in Figures 3- and 4 for various- values of 

t 

' b' . , From Fig. 3 it is verified that the, aperture distribution 
across the reflector is approximately Gaussian.' The- field at 
the rim is ; about 12 db below that in- the center. An example 
of , the transverse field distribution- is shown in Fig. 4. The 
mode, numbers were obtained by evaluation of the ..transverse 
field distributions. No correlation seems to exist between 
the mode, number n and the Q-value. The-.. only conclusion- that 
one can. draw at this- time- is . that the Q-value decreases- with 
increasing p. Further studies with modifications of the coupling 
elements between the resonator and the .waveguide- coupling energy 
in and out are in progress. 
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OPEN-WAVEGUIDE CHARACTERISTICS 
BY- FIELD- MATCHING- 

Abstract 

Field. matching method is, used, for considering a parallel- 
wall waveguide, Herewith the-, cross-section- is subdivided into 
regions,/ those 'between the. center conductor and the sidewalls 
and those 'adjacent to. the former. The- potential equations 
derived for the various regions contain two sets of unknown 
coefficients. To determine these coefficients, potential and. 
E-field are matched at the boundary surfaces between- the. re- 
gions* The- resultant set of infinite- simultaneous equations, 
is solved on the computer using a 20-term approximation. The 
potential ; distributions in x. and y directions along two lines 
are - plotted and compared with the result of a, numerical solution. 
The line impedance -is- obtained for several sizes of the-- center 
conductor and a plot is given. 



68 


Introduction , 

In this report, the field matching- method of solving elec- 
tromagnetic field an d^ wave propagation problems in parralel-wall' 
waveguides, is analyzed. It is studied in. relation to a rectangu- 
lar conductor- between the 'walls-. The 'results- are compared- with- 
the solution. of the static electric potential problem- using a 
simple- numerical technique. 

Field. Representation- 

Cons ider the parallel wall waveguide- with the inner con^- 

ductor shown in Figure- 1'. For TEM- operation 1 , Laplace 1 s._ equation 

is obtained- which- is applied in- the two regions I -‘and II: 

Y 

V & / = 0 -C^, d , CL 

V 4>Z = O C 

The two outer conductors 
are at zero potential and the 

Fig. 1 Cross-section 

inner conductor has normalized of waveguide 

potential of unity. 

Laplace 1 s . equation . is- 

solved in region- I by separation of variables under application 
of the boundary conditions: 

4>, -O @ &-o W a. ; as 

The solution may be- written as the. series - 





k A e* ALl 


d 
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In region II-, the solution to the potential is constructed 
by the following consideration: If the inner conductor is - ex- 

tended to infinity in both +x and -x directions , the potential 
in region II approaches the linear variation shown in figure 2 . 

<j> z = 2/c c ( y i c 

= / c ^ v< Vi 

This can be represented in 
a series of the form 

* ».!V 

in which the.B's are the Fourier Fig. 2' Potential distribution 

for infinite-width. 

coefficients- and may readily be- slab 

calculated. In the case of a slab of- finite- width, the potential 
distribution- in region II will deviate from the- above linear form 
by a function A<j> 2 . A<£ 2 is written as- superposition of two con- 

tributions, A^4> 2 and A 2 tf> 2 

A <f>z_ = A, + At 

A^ 4>2 is the deviation from the linear- distribution as- a result 
of removing the part of the. slab where x _< -d (Fig. 3-a) ; A 2 <j> 2 
is the. deviation- due to the removal of the part of the- slab where 
x >_ d' (Fig. 3-b) . 




Fig. 3 Cross section of guide with semi-infinite slab 





A^(J> 2 may be expressed by. a Fourier series 
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A i < P z - A? / 


_ inZT (V+cO 

e c % 


Similarly,. A 2 tf> 2 


is - expressed. as 


c „ J 












Note' that- A^cj >2 and A 2 <f> 2 are. equal to zero when- x approaches 

+ 00 - and -« respectively. Also, it is evident that- by symmetry, 

the. constants of the two series are- identically C The total 

m 

deviation is found to be 


A » A + AA 


where ' 


*£2y 

- . C*, Cos/S >nrrx mjt u 

, 

m — -Z (E . 


The total field in region II is 4> 2 = <f> 2 ' + A4> 2 which ma Y be - 
written in- either of the following- two forms r 


^ ^ ?<? + Cosh /L * '*£& 

n —ly 5 ^ v "* 




« s'# S c 


<?r 


tf- 


Or 


it 

At S - ^ 

h=t,3 r ~ 

1 


■+ 


c 
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Final Equations for Potential 
Rewriting the relations for the potential- in , each, region 
we- have; 

In region (1) , x >_ d 


$>, Cz-J) - 


dn <S. y 


° $ 9 v< £ 


In. region (2) , d> x >-d, the field may be written in- two forms 


(1) 


or 


^ fay) 

.X 


&Cx>yJ 

= 

S A? 


n=l,l~ 

< 3 . 

fay) 

= it 


4 fay) 

= i 



( 2 ) 


c- £ £ s< 


(3) 




Field Matching 

The coefficients A and, C are to be evaluated. This is 
■ n m 

done by matching the potential and the electric field at the 
boundary between the two. regions obtaining two systems of 
equations involving A'.s and C' s ; and- solving them simultaneously. 


Matching tjs: 

Evaluating and equating equations (1) and (2) or (1) and 
(3) at x = d we obtain 


2l a. a, 

A„ 



= ^ B„ AU »Hq 




C, <osh»*J Ait'S*', 

C e a 


(4) 
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or 




*«'>*,— 


M -<>*,' 


Z~ A„A^»JL« = / 

#=>,3,--- * ^ 


<? 


^ £ 3 v< *4 


(5 


nr 1 tt 


Each equation in (4) or (5) is ■ multiplied by sin y- and 

3. 

integrated in its proper, region-. 

Equation (4) gives: 

f C C 

J 2ZyAC*2£u Jy = J2_ C„ Cosh'nnd /u >Alu J<t 

<*/-L 

/ p . A " ' B " } ^ A “ d $ = o 

Adding the above and noting the orthogonality property a£ the 
left-hand side, we have: 


j « U -%) - 2 


^ P 


tyl ^ W*» 


(6 


where 


^ Cosh !2f<l J AC * « ACutHLy cJ% 


Equation (5) gives 


2_ + 


'• n ^>%— 


^4 
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Again adding the two equations and,, using the orthogonality on 
left-hand, side yields 


4 - 


a. 


l 


a. 


(»*) 


~h 




X 

He/A,'" 


c. 




( 7 ) 


which is equivalent to Eq. (6).- We "observe that the first right- 

hand term corresponds to. for the trapezoidal distribution 

according to Fig-. 2. Equations (6) or. (7) actually give only 

one system for a-, relationship between the A's and. C's. To find 

a second system of equations, another field quantity has to be 

matched at x = d. There . are two possibilities:, (a) matching 

E at x = d, (b) matching E at. x = d. 
x y. 


Matching E^: 

In each region 

In region (I) , 


4 


4 


3^ 





>.<? i % 


In region (II) 

.JL*? C su &*- ^ ^ 

Evaluating and. equating the two equations at x = d gives 




< C 


I 

Multiplying both, sides by sin —— y and integrating from 0 to c 
and' noting orthogonality property of the right-hand. side, "we - 
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obtain. 







^ SUZZJ . C„ 


where 




(3) 


Matching E^: 


In each' region E^. is-, found by: 


£u 


In region (I) 


2> 0 


£, -- t> 


~2J(x- cl) 


n7T 


4=43,~* **■ 


A„ e 


Coz*f 


£ 


In region (II) 


£,=±+^ 


***{*»-- 




^ £ *4 


r * <c 


= o 


«?{ V, 


Note that the result of obtaining the E y from (2) or (3) should 
be. equivalent. That is, the above set of equations- giving E^ 
in region II is equivalent to the relation obtained from (2): 


£t ^s?ic 


5 =27 ^ a s afi 1 


c < 7 $* % 
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Hence, in matching E at- x = ..d, either of the above representations 
may: be used. For completeness , let us. use both and find the re- 
sultant relations. 

Evaluate E y at x = d from both relations and. equate' the 


result. 






21 "S-?. £»!£!/ -t 


H ~ /;^ v • * 


°<U 


? Co,!, HI JCk ?<!■, 

^ c c c. 0 


( 9 ) 




or 


“ 2_ ?S.C« sefj, 

'X-. f A* i / oiyic 


h;/,X 


(10) 


51 If 4, fl* = 


< * < % 


Each equation in (9) is multiplied by - cos y and integrated 

cl 

in its proper region. 
c 


l ? « ■ A - 

^ 5 ? c ^ 0 ' 3 -V J f 

l ZLZ «.<*<*! e*#!'/ - * c . 




Add the two equations, use orthogonal properties in (0 - a) 
region- to obtain 


£ (a.-8,) 


where 


cr /r 

■= S '—fM / M 




an 




s' , hjrrsJ 
~2- Lo ^ h — 


Co$ ^ Cos ~Cy 
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n t m- 

If Eq. (10) -is multiplied by css — - y, , integrated in, proper 

ci 

region . and each part of it added., we obtain 



z Ck t > d 9 


f 






( 12 ) 


Summary. of the result of field matching is, as follows: 
Matching <j> at x = d: 


A (A„ - S„) = 


C P 

y, 


(Hz/jZj --- 


or 


4- <7 /!„ - ■ - — - AP nrr ^ ZT 

4 (»rr) z c c ~ 




Matching E^ at x = d: 


2 T 4* ■ G, 


^ t\u* 


M 

~2. 


I . c. 


Matching E^ at x = d: 


»r(A*-SL) - 21 C, • ^ 

tot!,*,-" 


/ 
tyt n 


or 


/- . / 

= 2_ <c '" 


( 6 ) 


(7) . 


( 8 ) 


(11) 


(12) 


To find A's and C's from the above, equation systems, two of the 
equations (6) , (7) , (8) , (11) , and (12) must be solved simultane- 
ously. Obviously Eqs . (6) and (7) taken simultaneously will not 
give any useful result-. Similarly (11) and (1-2) are not inde-- 
pendent. 
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The other possible- combinations which should be investigated are: 
Eq'S'. (6) or (7) - solved with Eq. (8), 

Eqs. (6) or (7) solved with Eqsv (11) or (12), 

Eqs-. (11) or- (12) solved with Eq. (8) . 


Solving [ (6) ,or (7)] wi-th f(H) or (12)]: 


4 - 

% 

- * >- 
^ tn 

C P 

'•— w HI n 

- 

% 

. ^ 2“ 

^ s 

2 : 

c m 

/ 2 /- 
L y? '/*?* 

J- /r ' ) 

<-7/7- 

Hr 




= 

MH 

Zosfi ££■< 6-/7 ^ »£? 



rr <*■ 




H1+! 

pL 

— 

mJT. C<?<,6 *zd. 

C- <L 

til 

rr a 


"V" v a * »rr A « c y 


= O 


to / c 


(»6) 3 --(»4x) 




(to/z) * - (»/#)*■ 



El O 


Therefore, no useful information is obtained. 


Solving [(6) -or (7)3 with C ( 8) [ r 

Find- from (8) and, substitute in. (7) 


c - 3- 


/ 


'“ 1 ~ torr <? ■ / hirrJ 

~c~~ yiz/ j5 ... 


$*, ' Q* 


* a A * “/TV— ^ — = E 

(krf) r « 


”Uli 

' * C 




/ S> . 'S- Qhtti > 






c 


(131 


i 3y-iy 
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This is solyed by taking a finite number of terms, then-A's- 
are - substituted- in (8) to find C's. 

Solving [ (8)] -with [(11) or (12)]; 

Again> C from Eg-. (8) is- substituted in Eq. (11): 
m -■ 




l 


2 

hirr 


t 

c. 








P 


if/ . 


or 


-5T A, ^ = -£ i7! % (1 

n -A/3 ," * ntz/,2,-. c ' 

k = /j 3 jP — 

Again., . the, above may- be expanded' by choosing a finite number 
of terms and the resulting system of equations, is solved for 
A's. The A's are .-substituted in either Eq. (8) or Eq. (11) 
to find C's. 

Computation of Capacitance and Impedance 
Solutions of.Eq. (13) were obtained for the cases where 
each series was, approximated by its first 20, 25, 30, 35 and, 40 
terms. It was observed that the value of the capacitance fluct- 
uated only in the 4th. significant figure (less than 0.1%). The 
results computed here are obtained from the, 2 0-term approximation 
of Eq. (13) in which case -indices k, n, and m are given the 
values 1,3,5 . . 39? 1, 3, 5 . . .39; and 1, 2, 3, . . .20 

respectively. 

Once the coefficients A's and C's are computed, the field 
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components may be determined by differentiating the potential 
functions. The total charge- per unit length of the guide is ob-- 
tained by integrating the normal component of the electric density 
vector at the . surface of the', guide walls; The capacitance is by 
definition the ratio of , total charge to the . potential difference.. 
Considering the symmetry of the structure, we write: 



Evaluating the integrals and. substituting. V •« 1, the capacitance 
is . 



The line impedance is 


where c 0 is the velocity of. light. - 

Fig. 4 shows, graph of line impedance versus c/a,, with 
2d/a as . parameter » 

Figures 5 and 6 show the comparison of the potential in the 
guide computed from the field-matching method and the numerical 
technique described- in the- appendix for one particular case where 
c/a=l<.l/4_, 2d/a = 1/2. It is observed that the two results agree 
quite well. The maximum deviations occur along the lines origin- 
ating from the comers of the inner conductor and in the- vicinity 


of the corner. 




Normalized potential 



Distance in y-directipn in a/40 units 


Fig. 5 Potential variations in y-direction , (the 

dots, indicate the result of field-matching 
method) . 


formalized Potential 



Distance in.. x-direction in a/40 units 

00 

Fig* 6 Potential variations in x-direction. (the 10 

dots, indicate the result of field-matching 
method ) d 
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APPENDIX 

• A numerical solution of the potential in the waveguide is . 

compared, with potential obtained from the field matching method,, 

The numerical technique outlined here is based- on finite-difference 

o 

solution of Laplace's equation. 

The potential near any point (x,y) may be expanded in a 
Taylor series disregarding higher-order terms 




A 

2- 

y- ^5: 
2_ 

~bx 2 

cfiCx-Ax^) 


-AclZ 

2- 

2. 

^0 

SIX 2 - 


4C*. ef+iy) S ?(*,}) y ^1? 
4c*. 9-4,) = 4 c*. y) -*9^ 


t 

-f 


-A?* 

2 


Combining the first two relations, we. obtain 


^ (x-Ax, % ) 



Ax 2 


Combining the next two equations gives 

ty+by-) ~ 2 + 4> (X, & - Ay) 

Ay* 




g 2 ^ gSi 

If — — and — =-2: are substituted in, the Laplace'^- equation and 


dx 4 


3y^ 


Ax = Ay = h, we -obtain 


= JL -b + 0(x ,%+*>) -h 0(x,if-h)J 

In order to utilize this relation, the cross-section of the 
waveguide is divided into a grid where the. size, of each mesh is 
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h x h-» According to the above equation, the potential at each 
node (grid intersection) . is obtained by averaging the -.-potentials 
of the four neighboring points (see. Fig„(7 ) * After the potentials 
of all ' the nodes which do not lie on boundaries - are ' successively 
computed, the process is repeated over again) each time using 
the previous values until each node-potential converges to a 
fixed value or. until its variations become, smaller than a pre- 
scribed tolerance c 

The- convergence- of the iterations may be improved by intro- 
ducing a procedure used in relaxation method for difference 
equations o The procedure. is to purposely under- or over-approximate 
the potential at each. step* If the change, of potential after the 
kth and (k + l)th- iteration is R 

k + 1 

then, instead; ©-£. finding <j> from 

,k+! r-,k 

<P = {ZT 4- £ 

we use 

A* ,k n ^k 

0 = <#> + ft £ 

where 8 may be less than or greater than unity and is. dependent 
on the geometry. 

As shown in , Fig 7, the- grid lines- are. chosen such that 
each boundary or. known potential surface coincides with a grid 
line. It should-, be noted that the accuracy of the converged 
potential values is dependent on the fineness of the. mesh. 

In Fig, 7 the boundary conditions are. given- by the zero 
potential on the walls at y = 0 and y = a. In addition, it is 
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assumed that the potential is- negligible- at +,x^ for x^ = 1.7 = 

This is justified since the potential- variation in the x-direction 
is proportional to exp{ — — — ) and for — > 1.5 its magnitude drops 
to less than- 0=01 of its value at the center conductor. 




[•1«] 


Black" Diagram of the Computer Program 
for the Numerical Solution 


Dimension array of 40 x 140 
P (40 , 140) 


.Set values of nodes on outer boundaries •= 0 


Set values of nodes on. inner conductor = ! 


Set values- of all'other nodes, by linear 
approximation from outer boundaries to 
inner conductor 


Starting at the first node, not- on the boundary, 
compute each node value by: 

R [P'(I-1,J)+P(I+1,J)+P(I,J-1)+P(I,J+1) ]-P(I,J) 

P (If J) » P(I, J) +" PR 

Skip the nodes, on the inner conductors,. 


Compute the value of next node 


All - nodes 
have not ,been 
computed 


All- nodes 
have been 
computed.- 


Check for 
convergence 
of values 


Not converged 


Converged 


Print out. the values of the 
nodes 
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